Finite-time collapse and soliton-like states in the dynamics of dissipative gases 
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A study of the gas dynamics of a dilute collection of the inelastically colliding hard spheres is 
presented. When diffusive processes are neglected the gas density blows up in a finite time. The 
blowup is the mathematical expression for one of the possible mechanisms for cluster formation in 
dissipative gases. The way diffusive processes smoothen the singularity has been studied. Exact 
localized soliton-type solutions of the gas dynamics when heat diffusion balances non-linear cooling 
are obtained. The presented results generalize previous findings for planar flows. 
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A simple case of a many-body dissipative system is pro- 
vided by a gas comprising particles whose collisions are 
inelastic. A common model considers a collection of hard 
spheres undergoing binary collisions, in which the normal 
component of the relative velocity of the particles is re- 
duced by a constant coefficient of restitution < r < 1 
(for r = 1 the elastic gas is recovered). This model of 
a granular gas is used to describe fluid-like motion of a 
dilute granular material - a collection of macroscopic par- 
ticles with inelastic interactions [J HI- Inelasticity models 
energy losses to the particles internal degrees of freedom. 
In this Letter we consider the dynamics of granular gases 
for r sufficiently close to 1, as explained below. 

At very small inelasticity (r < 1) the free behavior of 
a granular gas in a finite volume is very similar to the 
behavior of an ordinary gas. The gas evolution brings 
it to the spatially homogeneous state, differing from a 
molecular gas in the fact that the temperature decays in 
time due to inelastic energy losses. The decay law in this 
homogenous cooling state (below, the HCS) is a power 
with an inelastic cooling time-scale t c [3| • However as the 
inelasticity 1 — r is increased a "phase transition" (cf. 
occurs, still at 1 — r <C 1. The HCS becomes unstable 
and the gas starts developingclusters of particles (the 
clustering instability) ||J |f| [7|, [8| . 

Recently, under the condition l-r< 1 a derivation 
of the evolution towards the final state of a dilute granu- 
lar gas was given for particular vessel geometry [9| . Two 
assumptions were made: that the vessel is a channel and 
the channel length is not too large so that the sound 
travel time through the channel t s is much smaller than 
t c . In this geometry only longitudinal excitations of the 
gas arise which depend on the coordinate that measures 
distances in the long direction of the channel (the micro- 
scopic motions being in d = 2 or d = 3 where d is the spa- 
tial dimension) , while hydrodynamics is one-dimensional. 
The second assumption implies that sound travels back 
and forth along the channel many times before inelas- 
tic cooling becomes important. As a result, the pressure 
equilibrates on a time-scale of order t s while the pro- 
cesses associated with the inelasticity take place against 
the background of a spatially homogeneous pressure (cf. 
[H), which varies on the time-scale of the slow process t c . 
One can describe the slow dynamics using a single scalar 



equation. A study of this equation revealed that the non- 
linear development of the instability is very similar to the 
gas- liquid phase transition with many " droplets" arising 
throughout the gas which eventually coalesce to form a 
single drop of the new condensed phase. The final drop is 
described analytically as an exact soliton-type solution of 
the full system of gas dynamic equations. This solution 
has zero macroscopic velocity and a time-independent 
density profile. It describes a localized excitation of the 
gas. While the analogy between gra nula r ga s dynamics 
and phase ordering was discussed jldl . [Til [l2.Jl3j and ex- 
plored numerically [bj ] before, the theory in [0|is unique 
in that it provides a rather complete description of the 
evolution of the gas under the stated assumptions. In 
particular, it describes analytically the final state of the 
gas. 

The work [9| is the most recent in a series of papers 
d EE El, Cka H El devoted to the study of planar 
motions of a granular gas. The first study [l5[ intro- 
duced one-dimensional gas dynamics as a description 
of flows arising in a channel and studied the flows re- 
sulting from initial long scale excitations. The studies 
d, Il6t Il7l El, Gil were mainly devoted to ideal granu- 
lar gas dynamics (IGGD) obtained by neglecting heat 
conduction and viscosity. The rationale for studying 
the IGGD is that the clustering instability can be de- 
scribed in terms of the competition of two processes - the 
smoothing of inhomogeneities by viscous and thermal dif- 
fusion processes, having a characteristic time-scale tdif, 
and cooling dynamics corresponding to ideal, diffusion- 
less fluid dynamics, whose time-scale is t c . The HCS is 
unstable when the development of instabilities by cooling 
is faster than diffusion processes, t c < tdif. An interest- 
ing question is whether the density can grow indefinitely 
within the IGGD. It was found [17I. ll8L Il9j] that solutions 
of the one-dimensional IGGD become singular in a finite 
time t cr . The density field p diverges as p ~ l/(t cr — t) 2 . 
Remarkably, the gas pressure p = pT, where T is the 
temperature, remains finite at the singularity. The gas 
freezes at the singularity, T ~ {t cr — t) 2 , as can be roughly 
obtained from 9*T oc —pT 3 / 2 (see below) by substituting 
pT = const. Thus, at least for planar motions, the IGGD 
leads to an indefinite growth of the gas density. 

Below we show that the finite-time density singulari- 
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ties of the IGGD carry over to higher dimensions. The 
analysis is performed in the fast sound regime described 
above. Then we show that the soliton-type states exist 
also in higher dimensions and are not restricted to the 
fast sound regime. 

Consider both d = 2 (gas of discs in a plane) and d = 3 
(gas of hard spheres in a box) cases. Particles masses are 
set to unity and their diameter is designated by a. The 
IGGD equations read (see [13, [3 and references therein) 



d t p + (v ■ V) p = -pV • v, 
p [dtv + (v ■ V) v] = — Vp, 
9 ( T + (fV)T=-(7-l)TV 



ApT 



3/2 



(1) 

(2) 
(3) 



where v is the velocity, 7 is the adiabatic index of the gas 
(7 = 2 and 5/3 for d = 2 and d = 3, respectively), A — 
2 7r (<*-i)/2(l_ r 2) (T d-i/[ dr ( d /2)] (see e.g. [H) and T(. . . ) 
is the gamma- function. Equations C}-© differ from the 
hydrodynamic equations for a dilute gas of clastically col- 
liding spheres only by the presence of the inelastic cool- 
ing term —ApT 3 / 2 which is proportional to the product 
of the average energy loss per collision, ~ (1 — r 2 )T, and 
the collision rate, ~ pT 1 / 2 . The system ([I])-© applies 
when the conditions for the local thermal equilibrium are 
satisfied. Besides the usual conditions for gases [2l| , this 
leads to the quasi-elasticity condition 1 — r 2 <C 1. The 
inequality guarantees that the characteristic cooling time 
t c = 2/ Ap^ 2 p^ 2 inferred from Eq. ([3]) where po is the 
average gas density (the total gas mass divided by the 
volume of the channel), and po is a characteristic value 
of the initial pressure, is much larger than the mean free 



time T me an ~ VC ^ 1 Po Po )■ Then the system (P)-® 



describes well the evolution of long-scale excitations with 
large Reynolds number Re ~ pla d ~ 1 v/\/T, cf. [13] and 
Eqs. (fl"4|) - (fl"5"]) below. Here I is the excitation scale, p, v, 
T are characteristic values of the fields, and Re ^> 1 al- 
lows to neglect the heat conduction and viscosity terms 
in the equations. Below we use p instead of T as the 
second independent thermodynamic variable and employ 



d t p + (v ■ V)p = -7pV • v - Ap i/2 p 



l/2„3/2 



(4) 



instead of Eq. ((3|). We will assume small 1 — r so that 
the characteristic time-scale t s = L^/po/po, at which the 
sound traverses the characteristic vessel size L, obeys 
t s <C t c . Still we demand t s ^ T mean so that L is 
much larger than the mean free path l me an ~ l/poo"^ 1 - 
This leads to the range of vessel sizes l mean <C L <C 
lmean/ (1 — r 2 ). In this range the evolution of the gas is 
as follows: the gas relaxes to the state of local thermal 
equilibrium on time-scale T mean , then it relaxes to the 
state of a spatially homogeneous pressure on time-scale t s 
and only later, on timc-scalc t Cl the inelasticity becomes 
relevant. Thus the cooling dynamics occurs against the 
background of a spatially constant pressure, Vp = 0. In 
other words, the elimination of the fast acoustic modes 
allows so substitute the Eq. © by Vp = or p = p(t). 
To see this explicitly, we pass to rescaled variables. We 



measure time in the units of t c , coordinate in the units of 
L, and the gas density, pressure and velocity in the units 
of po, Po and L/t c respectively. Keeping the original no- 
tation for the rescaled variables, we observe that Eq. ([T|) 
does not change, while Eqs. ([2]) and (|4]) become 



(t s /t c ) 2 p [d t v + (v ■ V) v] = -Vp, 
d t p+ (v-V)p = -7pV • v - 2p 1/2 p 3/2 - 



(5) 
(6) 



The leading order approximation in t s /t c < 1 is ob- 
tained by dropping the velocity term in the Euler equa- 
tion which results in Vp = 0. Thus the pressure depends 
on time only, p = p(t). Averaging Eq. ([B]) over the volume 
and using v ■ Vp = we find 



P (t) = -2(pV 2 (t))pV 2 (t), 



(7) 



where angular brackets stand for the spatial average, 
(p 1 / 2 ^)) = / p l t 2 (x,t)dx/V with V the system volume. 
We used that (V-u) =0 both for periodic or zero bound- 
ary conditions on velocity. Expressing now V • v from 
Eq. ^ in terms of p and p, we rewrite Eq. |T]) as 



livpr 1 ' 2 [d t +v-V]p 



1/2 



P 



1/2 



(P 



l/2\ 



(8) 



The above reduction of the original system of Eqs. ([T])- 
(@| has great simplicity: the density field is driven by 
the spatial fluctuations of p x / 2 {x,t). It is possible to 
write down a closed equation for the average (p 1 ^ 2 )- We 
note that Eq. {J) can be written as dtp 1 ! 2 + (v ■ V)p^ 2 = 
— p 1 / 2 V-t>/2. Averaging the last equation over space and 
using that integration by parts implies — (p 1 / 2 V ■ v) = 
((u-V)p 1 / 2 ), we find ((v-V)p 1 / 2 ) = -2d t (p 1/2 }. Finally, 
performing spatial averaging of Eq. ([8]) multiplied by p 1 / 2 
and using (p) = 1 we obtain 



d X /dr = -l + x 2 , X (t) = (p 1/2 )(t), 



(9) 



where following Q we introduced a new time variable 

t = J *p 1 / 2 (t')<it'/7- The last two equations, though 
they cannot be derived by the same means as in the one- 
dimensional case, have exactly the same form as in that 
case [9( . The solution of Eq. with the initial condition 
Xo = X(t = 0) is 



X(r) 



Xo 



tanhr 



1 — Xo tanh ' 



(10) 



The Cauchy-Schwarz inequality (fg) 2 < (f 2 )(g 2 ) applied 
to / = p 1 / 2 and g = 1 gives (p 1 ' 2 ) < (p) 1 ^ 2 - Combined 
with (p) = 1 this implies that xo ^ 1 with equality hold- 
ing only in the case of spatially homogeneous density. As 
a result the above expression for x( T ) becomes negative 
at some finite time tq. Since by definition x( T ) > this 
signifies that the solution must break down at r < tq. 
Indeed the density evolving according to Eq. ([S]) blows 
up at t < tq. To see the blowup we note that the den- 
sity p(r, a) in the (Lagrangian) frame moving with the 



3 



fluid p(t,o) = p[r, x(t, a)] with d T x(r,a) = v(t,x(t,o) 
obeys (we omit tildes with no ambiguity) 

p-^ir, a)d T p^ 2 (r, a) = p^ 2 (r, a) - X (r), (11) 

where a is a label of the trajectory (say the initial coor- 
dinate x(0, a) — a). Following [§], we pass to w = l/^/p, 
that obeys a linear equation 



O t w(t, a) — x{t)w(t, a) 



-1. 



(12) 



The above equations are the same as in one-dimensional 
situation and one can take the solution from Q. The 
density in the Lagrangian frame is given by 



w(t, a) 



w(0, a) + xo [cosh r — 1] — sinh r 
cosh t — xo sinh r 



(13) 



The solution describes the formation of density singular- 
ity exactly as in d = 1 situation, see details in [§]. The 
only difference is that in d > 1 the knowledge of den- 
sity in the Lagrangian frame does not allow to recover its 
profile in real space, while in d = 1 this is possible Q. 
Thus we know that the singularity occurs and we know 
its law, but we don't know where in space it occurs. The 
resolution of this question demands the use of vorticity 
equation for the calculation of v{x,t), see [23j], and it is 
beyond the scope of the present work. 

Thus we have shown that the evolution of any inhomo- 
geneous initial condition according to Eqs. CO)-® pro- 
duces finite time singularity of the density for vessel sizes 
Imean <C L -C l m ean/(^ ~ i" 2 ) . We expect the finite time 
blowup holds in the whole range of applicability of hy- 
drodynamics Imean <C L. The reason is that the system 
of Eqs. (HJ-© is hyperbolic so that creation of singular- 
ities is determined by spatially localized dynamics that 
fits the previous analysis. This viewpoint is also con- 
firmed by the existence of a group of exact solutions for 
planar flows that have the same finite-time singularity 
|l~7l [H[ . Still further work is needed to establish the fact 
completely. 

Does the heat conduction and viscosity arrest the 
above singularity? The answer is not obvious due to the 
vanishing of heat conductivity and viscosity at a singu- 
larity with T — 0. Here we show indirectly that the con- 
duction arrests the singularity by displaying exact, valid 
at any L, solutions of equations with heat conduction 
for which the density is not growing. Diffusive processes 
add to Eqs. (J2J) and ^ the terms describing the viscous 
momentum transfer and the conduction of heat: 



p [d t Vi + (v ■ V) Vi] = -Vip 

d t p + (v ■ V) p = — 7pV • v - 
+ (2k /3)V 2 T 3 / 2 + ( 7 - l^oVTaijdjVi, 



Ap^p 3 / 2 



■(14) 



(15) 



diVj - 



while Eq. |T]) is unchanged. Above cry = djVi 
2SijV ■ v/d and it is enough for our purposes here that 
the coefficients v$ and kq are of order l/cr d_1 Let 
us search for the solutions of the above system for which 



the gas performs no macroscopic motion, v = 0. Then 
Eqs. H]) and {Tj) produce p{x,t) = p(x) and p(x,t) — 
p(t). Further, spatial averaging of Eq. (fT5|) gives p = 
— Ap 3 / 2 (p 1 / 2 ) which substitution back in Eq. (fT5]) gives 
the equation defining the stationary profile p(x): 



p l/2 _ {p l/ 2) _ (1/3)V 2 /0 -3/2 = Q _ 



(16) 



where we measure p in the units of po and distances 
in the units of l cr = \/2kq / Ap^. The scale l cr is the 
scale at which the cooling and the heat diffusion bal- 
ance each other and it equals the critical length be- 
yond which spatial perturbations of the HCS are unsta- 
ble 0, 0]. The demand of zero heat flux through the 
boundary and T(x,t) = p(t)/ p(x,t) lead to the Neu- 
mann boundary condition (vanishing normal component 
of Vp) for Eq. (fT6|) . The equation has inhomogeneous 
solutions which physics is the same as in other soliton 
solutions of non- linear physics (see e. g. [Hj]): diffu- 
sion and non-linearity balance each other. The solutions 
for d = 1 were considered in Q. It was shown there 
that inhomogeneous solutions of Eq. (|16|) exist only for 
L > l cr (in contrast to the solution p = const that exists 
always). As a result at L > l cr where the HCS is unstable 
the soliton-type solution becomes the global attractor of 
the dynamics 0. The solutions however exist also for 
d > 1. As an example consider the d — 3 spherically 
symmetric solution of Eq. (| 1 6|) with p = p(r) obeying 
p'(0) = p'(L) = 0. The numerical solution of the problem 
is shown in the Figure. As in d = 1 case inhomogeneous 




FIG. 1: Numerical plot of $ e 
The system length L obeys L 



I l/2\3 -3/2 

(p ' ) p ' versus x 

— 6/cr- 



r/L. 



solutions exist only for L > l cr (it can be shown that l cr 
corresponds to the same length for which the homoge- 
neous cooling state becomes unstable versus spherically 
symmetric perturbations). Note that p = pT is propor- 
tional to the energy density since pv 2 = 0. One has 
p = po[l + A(p 1 / 2 }p^ 2 t/2]~ 2 both for homogeneous and 
inhomogeneous density solutions however for inhomoge- 
neous solutions the decay is slower by {p 1 ^ 2 ) < (p) 1 ^ 2 , 
cf. [!|. Thus the transition from the HCS to the soli- 
ton state is the transition to a state with a slower energy 
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decay rate. The above solutions are general and do not 
involve the fast sound assumption. 

It is natural to ask if the soliton solutions are global 
attractors as in d — 1 case. One can start with the case 
where sound is the fastest macroscopic process in the sys- 
tem. Then again the evolution occurs against the back- 
ground of a spatially uniform pressure. One finds that 
thermal conduction effects lead to an additional, Lapla- 
cian, term in the effective equation ([8j): 

Clearly soliton solutions are also stationary solutions of 
the above effective dynamics. In d = 1 introducing the 
Lagrangian mass coordinate m(x,t) = p(x',t)dx' one 
may get a closed description for the evolution of w(m, t). 
Eq. (JT2J) gets modified to (note that m is a particular 
case of a Lagrangian label a) 

d T w(T,m) — x{t)w(t,itl) = —l + ur 1 (T,m)df n w(T,m). 

It is the study of the above equation that has allowed the 
complete understanding of the gas evolution in d = 1 • 
In d > 1 such reduction is impossible (derivatives of La- 
grangian coordinates with respect to the initial position 
are no longer expressible via density only). The study of 
the above dynamics is a subject for future research. 

The above results generalize previous findings for pla- 
nar hydrodynamic flows [jj [l?], [H| showing that the anal- 
ysis in d = 1 can serve as a good guide in the study of the 
fully three (or two) - dimensional hydrodynamics of gran- 
ular gases. The non-linear ideal dynamics of free granular 
gases, obtained by neglecting heat conduction and viscos- 
ity, leads to density divergence in a finite time just like in 



d = 1. This singularity provides a robust distinction be- 
tween the dissipative and conservative gases (where the 
density always remains bounded) and describes a possi- 
ble mechanism for the cluster formation. The demon- 
stration was performed in a particular range of vessel 
sizes and it was argued that the hyperbolic nature of the 
equations implies that the blowup holds generally. The 
property that the arising singularities are effectively one- 
dimensional is known for a simpler gas dynamics where 
the velocity obeys the Hopf equation d t v + (v ■ V)t> = 
and the density obeys the continuity equation. In both 
cases, the sign of effective one-dimensionality of singular- 
ities is the arising filament structure of the high density 
regions - the singularities occur in the direction perpen- 
dicular to filaments. The filament structure of high den- 
sity regions in granular gases is observed in [5|]. 

We derived new soliton-like states of the gas as exact 
solutions of the full system of gas dynamic equations. 
The solutions show how heat conduction smooths the 
density singularity described above. They generalize d = 
1 solutions that were shown to describe the final state of 
the gas at least in the fast sound regime 0]. Are soliton 
solutions stable and do they describe the final state of 
granular gas in d = 2 and d = 3? Does gas evolution 
parallel the gas- liquid transition as in d = 1? These are 
subjects for further work. 
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